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We investigate the problem of determining the parameters that describe a quantum channel. 
It is assumed that the users of the channel have at best only partial knowledge of it and make 
use of a finite amount of resources to estimate it. We discuss simple protocols for the estimation 
of the parameters of several classes of channels that are studied in the current literature. We 
define two difi'erent quantitative measures of the quality of the estimation schemes, one based on 
the standard deviation, the other one on the fidelity. The possibility of protocols that employ 
entangled particles is also considered. It turns out that the use of entangled particles as a new 
kind of nonclassical resource enhances the estimation quality of some classes of quantum channel. 
Further, the investigated methods allow us to extend them to higher dimensional quantum systems. 

o 
o 

CsJ ■ I. INTRODUCTION 

^ ' Quantum information processing has attracted a lot of interest in recent years, following Deutsch's investigations 
' 101 concerning the potentiality of a quantum computer, i.e., a computer where information is stored and processed in 
quantum systems. Their application as quantum information carriers gives rise to outstanding possibilities, like secret 
communication (quantum cryptography) and the implementation of quantum networks and quantum algorithms that 

^ ^ are more efficient than classical ones |^ . 

^ ' Many investigations concern the transmission of quantum information from one party (usually called Alice) to 
f — , another (Bob) through a communication channel. In the most basic configuration the information is encoded in 
CO ' qubits. If the qubits are perfectly protected from environmental influence. Bob receives them in the same state 
', prepared by Alice. In the more realistic case, however, the qubits have a nontrivial dynamics during the transmission 
because of their interaction with the environment Q|. Therefore, Bob receives a set of distorted qubits because of the 
disturbing action of the channel. 

Up to now investigations have focused mainly on two subjects: Determination of the channel capacity Q and 
reconstruction schemes for the original quantum state under the assumption that the action of the quantum channel 
is known Here we focus our attention on the problem that precedes, both from a logical and a practical point of 
view, all those schemes: The problem of determining the properties of the quantum channel. This problem has not 
been investigated so far, with the exception of very recent articles |^,^. The reliable transfer of quantum information 
■ requires a well known intermediate device. The knowledge of the behaviour of a channel is also essential to construct 
^ quantum codes 0J^. In particular, we consider the case when Alice and Bob use a finite amount N of qubits, as this 
is the realistic case. We assume that Alice and Bob have, if ever, only a partial knowledge of the properties of the 
^ ' quantum channel and they want to estimate the parameters that characterize it. 
• i-H , The article is organized as follows. In section ^ we shall give the basic idea of quantum channel estimation 
' and introduce the notation as well as the tools to quantify the quality of channel estimation protocols. We shall 
^ J then continue with the problem of parametrizing quantum channels appropriately in section [fn|. Then we are in a 
- - ■ position to envisage the estimation protocol for the case of one parameter channels in section PV|. In particular, we 
shall investigate the optimal estimation protocols for the depolarizing channel, the phase damping channel and the 
amplitude damping channel. We shall also give the estimation scheme for an arbitrary qubit channel. In section ^ 
we explore the use of entanglement as a powerful nonclassical resource in the context of quantum channel estimation. 
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section VI deals with higher dimensional quantum channels before we conclude in section VII 



II. A GENERAL DESCRIPTION OF CHANNEL ESTIMATION 



The determination of all properties of a quantum channel is of considerable importance for any quantum commu- 
nication protocol. In practice such a quantum channel can be a transmission line, the storage for a quantum system, 
or an uncontrolled time evolution of the underlying quantum system. The behaviour of such channels is generally not 
known from the beginning, so we have to find methods to gain this knowledge. 

This is in an exact way only possible if one has infinite resources, which means an infinite amount of well prepared 
quantum systems. The influence of the channel on each member of such an ensemble can then be studied, i.e., the 
corresponding statistics allows us to characterize the channel. In a pratical application, however, such a condition 
will never be fulfilled. Instead we have to come along with low numbers of available quantum systems. We therefore 
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cannot determine the action of a quantum channel perfectly, but only up to some accuracy. We therefore speak of 
channel estimation rather than channel determination, which would be the case for infinite resources. 

A quantum channel describes the evolution affecting the state of a quantum system. It can describe effects like 
decoherence or interaction with the environment as well as controlled or uncontrolled time evolution occuring during 
storage or transmission. In mathematical terms a quantum channel is a completely positive linear map C (CP-map) 
|^,|o|, which transforms a density operator p to another density operator 

p' ^Cp. 

Each quantum channel C can be parametrized by a vector A with L components. For a specific channel we shall 
therefore write throughout the paper. Depending on the initial knowledge about the channel, the number of 

parameters differs. The goal of channel estimation is to specify the parameter vector A. 

The protocol Alice and Bob have to follow in order to estimate the properties of a quantum channel is depicted 
in figure |l|. Alice and Bob agree on a set of N quantum states pi,i = 1,2,. ..,A^, which are prepared by Alice and 
then sent through the quantum channel C^. Therefore, Bob receives the N states p'^ = pi. He can now perform 

measurements on them. From the results he has to deduce an estimated vector A"*'* which should be as close as 
possible to the underlying parameter vector A of the quantum channel. 




J 

FIG. 1. Basic scheme of channel estimation. Alice sends A*' quantum state pi to Bob. The channel maps these states onto 
the states p'^ — C^pi, on which Bob can perform arbitrary measurements. Note that Bob's measurements are designed with the 
knowledge of the original quantum states pi. His final aim will be to present an estimated vector A"''' being as close as possible 
to the underlying parameter vector A. 



How can we quantify Bob's estimation? To answer this we introduce two errors or cost functions which describe 
how good the channel is estimated. 

The first obvious cost function is the statistical error 

c.{N,X)^J2i^^-^THN)f (1) 

in the estimation of the parameter vector A. Note that the elements of the estimated parameter vector A''*'* strongly 
depend on the available resources, i.e. the number N of systems prepared by Alice. We also emphasize that Cs 
describes the error for one single run of an estimation protocol. However, we are not interested in the single run error 
but in the average error of a given protocol. Therefore, we sum over all possible measurement outcomes J to get 
the mean statistical error 

MN,\) = {cs{N,X))j 

while keeping the number N of resources fixed. 

Though this looks as a good benchmark to quantify the quality of an estimation protocol it has a major drawback. 
The cost function strongly depends on the parametrization of the quantum channel. While this is not so important 
if one compares different protocols using the same parametrization it anyhow would be much better if we could give a 
cost function which is independent of any specifications. We define such a cost function with the help of the average 
overlap 
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T{C,,C2)^ (F{C,\4,){^\,C2\4^)m) (2) 



between two quantum channels C\ and C2, where we average the fidehty flT 



F(pi,p2) = Tr^^VprP2\/pr (3) 

between two mixed states pi ^ Ci p and p2 = C2P over all possible pure quantum states p = |V')('0l Q- Since the 
fidelity ranges from zero to one, the fidelity error is given by 

Cf{NS)^l-HCl.Cy^...^M)) (4) 

which now is zero for identical quantum channels. Again we average over all possible measurement outcomes to get 
the mean fidelity error 

Cf{N,X)^{cf{N,X))j (5) 

which quantifies the whole protocol and not a specific single run. 

In the first part of this paper we are only dealing with qubits described by the density operator 

p=^il + s-d) (6) 

with Bloch vector s — Tr(p(T) and the Pauli matrices a = {(JxT<^yT<^z)- The action of a channel is then completely 
described by the function 

s'^C{s), (7) 
which maps the Bloch vector s to a new Bloch vector s' . In particular, the mixed state fidelity equation (|3|), for two 



qubits with Bloch vectors si and 5*2, see equation (|6|), then simplifies to |11 



, = - ^1 + Si . S2 + ^{l-\s,\^){l-\s2V)\ (8) 
which leads to an average channel overlap, equation (|^), 

HCuC^)^ J ^ F(Ci(7i),C2(n)) (9) 

for the two qubit channels Ci and C2. As emphasized above we only average over pure input states with unit Bloch 
vector 

fi = (cos $ sin 0, sin $ sin 8, cos 8) (10) 

and Bloch sphere element dVl = sin 8 dQ d^. By inserting equations (^) and (^ into equations (^) and (^) we get a 
cost function for the comparison of two qubit channels which is independent of the chosen parametrization. 

III. PARAMETRIZATION OF A QUANTUM CHANNEL 

As we have already mentioned in section ||, the qubits that Alice sends to Bob are fully characterized by their 
Bloch vector. Therefore, the disturbing action of the channel modifies the initial Bloch vector s according to equation 
(^. It has been shown that for completely positive operators the action of the quantum channel C is given by an 
affine transformation 

g = Ms + v (11) 



^ Since for our estimation protocol we are only sending pure states through the quantum channel we also only average over 
pure states. 
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where M denotes a 3 x 3 invertible matrix and w is a vector. The transformation is thus described by 12 parameters, 
9 for the matrix Ai and 3 for the vector v. These 12 parameters have to fullfill some constraints to guarantee the 



complete positivity of C |12|. The definition of the parameters is somewhat arbitrary. We wiU not use the parameters 
as defined in |12|, but adopt a different parametrization 
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(12) 



in terms of the parameter vector A — (Ai, . . . , Ai2)"^. In general, the choice of the best parametrization of quantum 
channels depends on the relevant features of the channel and also on which observables can be measured. The reason 
for the choice of the parametrization, equation ([ij) , will become clear at the end of the section [V , in which a protocol 
for the general channel is described. 



IV. ESTIMATION OF CHANNEL PARAMETERS 



Although the characterization of the general quantum channel requires the determination of 12 parameters, only 
a smaller number of parameters must be determined in practice for a given class of quantum channels. Indeed, 
the knowledge of the properties of the physical devices used for quantum communication gives information on some 
parameters and allows to reduce the number of parameters to be estimated. We shall now examine in detail some 
known channels described by only one parameter. 



A. The depolarizing channel 

The first channel we consider is the depolarizing channel. The relation (|l^) between the Bloch vectors s and s ' of 
Alice's and Bob's qubit, respectively, reduces to the simple form 

SA ' = (1 - 2A)s, 

where < A < 1/2 is the only parameter that describes this quantum channel. This channel is a good model when 
quantum information is encoded in the photon polarization that can change along the transmission fiber via random 
rotations of the polarization direction. If Alice prepares the qubit in the pure state p — \'>p){'ip\, the qubit Bob receives 
is described by the state 

p' = (l-A)|V)(^|+A|V;)(V5|, (13) 

where j-i/i) denotes the state orthogonal to lip) . We note that the depolarizing channel has no preferred basis. Therefore, 
its action is isotropic in the direction of the input state. This means that the action of the channel is described by 
equation (^3|) even after changing the basis of states. 

Bob must estimate A, which ranges between (noiseless channel) and 1/2 (total depolarization). The estimation 
protocol is the following: Alice prepares N qubits in the pure state 

|k)-e-*/2cos||U+e'*/2sin||U 

which denotes spin up in the direction n, equation ([l0|), in terms of eigenvectors |tz) and ||z) of (Tz- The state |tn) 
is sent to Bob through the channel. Bob knows the direction n and measures the spin of the qubit he receives along 
ft. The outcome probabilities are 

P(Tfl) = l-A, P(k) = A. 

Since Alice sends a finite number of qubits, Bob can only determine frequencies of measurements instead of probabil- 
ities. After Bob has measured i < N qubits with spin down and the remaining N — i with spin up, his estimate of A 
is X'^^ = i/N. Note that in a single run of the probabilistic estimation method we can get results X°^^ > 1/2, which 
is outside the range of the depolarizing channel. However, for the calculation of the average errors we do not have to 
take that into account. The cost function Cg, equation (11^), is thus 
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The average error is easily obtained when one considers that Bob's frequencies of measurements occur according to a 
binomial probability distribution, since each of the i measurements of spin down occurs with probability A and each 
of the spin up measurements occurs with probability 1 — A. Therefore, the mean statistical error reads 



N 



[N, A) = ^ (^)a'(I-A)^-^ (A-A^^*)' 



N 



i=0 



N 



N 



A(l-A) 
N ' 



(14) 



This function, shown in figure ^ a), scales with the available finite resources N. It vanishes when the channel faithfully 
preserves the polarization, A = 0. The largest average error occurs for A = 1/2, when the two actions of preserving 
and changing the polarization have the same probability to occur. 
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FIG. 2. a) The average statistical mean error Cs{N,X) equation ([l4|), and b) the fidelity mean error Cf(N,\), equation (p^, 
for a depolarizing channel with parameter A, are shown for different values of the number of qubits N. 



The cost function c/ is obtained from equations and (^) 



c,{N,\) = 1- 



which leads to the fidelity mean error 

N 



xf~'cf{N,X) 



c-,(7V,A)=^(^)a^(1 



(15) 



which is shown in figure |2|b). Although c/ does not show a simple 1/A^-dependency, it clearly decreases for increasing 
values of the number of qubits N. For a given quality Cs or c/ one can use figure || and read off the number iV of 
needed resources. 



B. The phase damping channel 



The phase-damping channel acts only on two components of the Bloch vector, leaving the third one unchanged: 
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sa' = ((1-2A) si,(l-2A) S2,s3f . 

Here < A < 1/2 is the damping parameter. This channel, contrarily to the depolarizing channel, has a preferred 
basis. In terms of density matrices, it transforms the initial state 



where = (l^l p |tz), etc., into 



_ I Pll Pit 

Pn Pn 



Pll (l-2A)pxt 



(16) 



We note that here the parameter A only appears in the off-diagonal terms. For this reason, the phase damping channel 
is a good model to describe decoherence |lj] . Indeed, a repeated application of this channel leads to a vanishing of 
the off-diagonal terms in p' , whereas its diagonal terms are preserved. 

Since the parameter A of the phase damping channel appears in the off-diagonal elements of p' , the protocol is 
the following: Alice sends N qubits with the Bloch vector in the x — y plane; for instance, qubits in the state 
ITi) = (Uz)+ |Tz))/V2 can be used. This would correspond to a density operator whose matrix elements are all equal 
to 1/2, can be used. The density matrix of Bob's qubit is then given by 



1 (1 - 2A) 
(1 - 2A) 1 



Now he measures the spin in the x direction. The theoretical probabilities are 

P{U) = {U\ p' = I - \ 
P{U) = iU p' \U) = X, 

and we denote the frequency of a spin-down measurement as i/N, leading to A*"** = i/N. The mean statistical error 
has again the form 



TV 



i=0 



c-.(7V,A) = 5]('y)A'(l-A) 



A(l - A) 
N ' 



as for the depolarizing channel. For the fidelity, equation (||), we obtain 

-F(sa',5 lost) - 1 + [2AA'=^' - A - A'^'^' 



+2v/A(l-A)A<'«t(l- A°^') (sl + sl) 



(17) 



and after averaging over the Bloch surface, according to equation 



Therefore the fidelity mean error c/ reads 



1 - 3A - -A'^^' + -AA^^* + -VA(1-A)A-H1-A-') 



N 



A+(1-2A)--2JA(1-A)- 1-- 



N 



N 



N 



Ml-A)--5]f )y{l-Xf~WKl-mN-i) 



(18) 



which is shown in figure y b). This mean error is very similar to that obtained for the depolarizing channel. 
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FIG. 3. a) The average statistical mean error Ca{N, A) , equation (|l7|), and b) the fidelity mean error Cf{N, A), equation (|l8|), 
for a phase damping channel with parameter A, are shown for different values of the number of qubits N. 



C. The amplitude damping channel 



The amplitude damping channel affects all components of the Bloch vector according to 



sx ' = (yi - A si, Vl - A S2, (1 - A) S3 + A 
where < A < 1 is the damping parameter. The density matrix, equation (n 



(19) 



is transformed into 



A+ (1- A)pa VI- Apxt 
VI - Apix (1-A)/7|| 

This channel is a good model for spontaneous decay from an atomic excited state \1z) to the ground state \lz)- 
Repeated applications of this channel cause all elements but one of the density matrix to vanish. Now the parameter 
A appears in all the elements of p' and the channel clearly possesses a preferred basis. 

If Alice and Bob know that they are using an amplitude damping channel, Alice sends all N qubits in the state 
Itz). The density operator of the qubit received by Bob is 

A 

1-A 

He measures the spin along the z direction (the diagonal elements of p' are the probabilities to find spin down and 
spin up, respectively). We denote the frequency of spin down measurements with i/N. The statistical cost function 
is again 



N 



A 



N 



A(l-A) 

TV 



(20) 



Using equation (19) we obtain 



F{sx ', s ',eeO = - [l + V(l-A)(l-A-') (si + si) + (1 - A)(l - A-*)s2 

+ [A(1-A^^*) + A^^'(1-A)] S3 
+AA<='^* + (1 - S3)VA(1- A)VA'=^H1- A°^') 



for the fidelity cost function, equation (S), and 



^(CA,CAesO = - [4 + 2v/(l-A)(l-A-')-A-A^ 
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for the averaged fidelity (m. Thus 
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(21) 



which is illustrated by figure ^. 
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FIG. 4. a) The average statistical mean error Ca{N, A) , equation (^), and b) the fidelity mean error c/ (iV, A), equation (|l|), 
for the amplitude damping channel with parameter A, are shown for different values of the number of qubits N. 

We see from the figures (||)-(||) that for a fixed number of resources the mean error c/ has a local maximum for 
small values of the parameter A. For the amplitude damping channel there is also a second local maximum for large 
values of A. 

D. The general quantum channel 

We now come to the problem of estimating the 12 parameters A = (Ai, . . . , Xi2)^ of the general quantum channel. 
The protocol is summarized in table 1: Alice prepares 12 sets of qubits, divided into 4 groups. Bob measures the 
spin of the three sets in each group along x, y, and z, respectively. From the measurements on each set he gets an 
estimate of one of the parameters A^. The parametrization (|l^ has been chosen for this purpose. 

The statistical cost function for the general channel is thus a generalization of the cost functions for one-parameter 
channels. 



N/12 N/12 

il=0 412=0 

12 

= ^— VA,(1- A,,) 

' .7 = 1 



n—1 ^ ^ 



The mean fidelity error can be calculated numerically. However, we do not give the expression here as it gives no 
particular further insight. 
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TABLE I. The protocol for the general channel. Alice prepares 4 sets of qubits each in a state corresponding to the 
spin pointing in a well prepared direction (left); Bob splits each set into 3 subsets of N/12 qubits each and performs a spin 
measurement along x, y or z, respectively (center); the frequency of each set of spin measurements gives the estimate of one 
channel parameter (right). 
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V. THE PAULI CHANNEL FOR QUBITS 



Up to now we have only considered estimation methods based on measuring single qubits sent through the quantum 
channel. However, we are not restricted to these estimation schemes. Instead of sending single qubits we could use 
entangled qubit pairs In this section we will demonstrate the superiority of such entanglement-based estima- 
tion methods for the estimation of the so called Pauli channel by comparing both schemes (for a different use of 
entanglement as a powerful resource when using Pauli channels, sec p^). 

The Pauli channel is widely discussed in the literature, especially in the context of quantum error correction [Q. 
Its name originates from the error operators of the channel. These error operators are the Pauli spin matrices ax^ (Jy 
and (Jz- The operators define the quantum mechanical analogue to bit errors in a classical communication channel 
since Ux causes a bit flip (| |) transforms into 1 1) and vice versa), cr^ causes a phase flip (| i) + | T) transforms into 
I I) — I t)) and ay results in a combined bit and phase flip. The Pauli channel is described by three probabilities 
A = (Ai, A2, Aa)"^ for the occurrence of the three errors. If an initial quantum state p is sent through the channel, the 
Pauli channel transforms p into 

p' = (1 - Ai - A2 - \s.)p + \i<JxP<Jx + \2Oypoy -f Aacr^pcTz. (22) 

Thus we see that the density operator p remains unchanged with probability 1 — Ai — A2 — A3, whereas with probability 
Ai the qubit undergoes a bit flip, with probability A3 there occurs a phase flip, and with probability A2 both a bit flip 
and a phase flip take place. 



A. Estimation using single qubits 



We will flrst describe the single qubit estimation scheme for the Pauli channel, before we switch to the entanglement- 
based one in the next subsection. Following the general estimation scheme presented in section || the protocol to 
estimate the parameters of the Pauli channel, equation (p2|), requires the preparation of three different quantum states 
with spin along three orthogonal directions. Alice sends(i) M = iV/3 qubits in the state | "f^), (ii) M qubits in the 
state I ta;), and (iii) M qubits in the state | "fj,). Bob measures their spins along the direction of spin-preparation, 
namely the z, x, and y axes, respectively. The measurement probabilities of spin down along those three directions 
are then given by 

P(i.) = A1 + A2, 
P(i.) = A2 + A3, 
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P{iy) = M + X3 



The estimated parameter values can be calculated via 
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where ik/M, k = 1,2,3, denote the frequencies of spin down results along the directions x, y and z. We note that, 
although the probabilities are positive or vanish, their estimated values A^*^* may be negative. This occurs because 
in the present case the measured frequencies are not the estimates of the parameters. Nonetheless, the average cost 
functions can always be evaluated. For the statistical cost function we find 
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= ^ [^1(1 - ^1 - ^2) + A2(l - A2 - A3) 
+A3(l-A3-Ai)]. 



(23) 



for the average error c/°p of the estimation with separable qubits |6) . For fixed N the average error has a maximum 
at Ai = A2 = A3 = 1/4, in which case all acting operators occur with the same probability. On the other hand, the 
average error vanishes when faithful transmission or one of the errors occurs with certainty. 

Instead of using the statistical cost function Cg we can also use the fidelity based cost function c/ (jj). The average 



error c 



of the estimation via separable qubits is then given by 

Ar/3 N/3 N/3 



— sen 



(iV,A) 



Ei:i:(T)(T)(T)'^'«'"'-^"'»""- 

n=0 (2=0 13=0 ^ i/V / \ 3 / 

xiPiuyrH^-piu))''^'-'' 



(24) 



The cost function c/ cannot be calculated analytically. We shall compare the results (^3|) and ( |24[ ) with the same 
mean errors for a different estimation scheme, where entangled pairs are used, that we are now going to illustrate. 



B. Estimation using entangled states 



All the protocols that we have considered so far are based on single qubits prepared in pure states. These qubits 
are sent through the channel one after another. However, one can envisage estimation schemes with different features. 
An interesting and powerful alternative scheme is based on the use of entangled states . In this case the estimation 
scheme requires Alice and Bob to share entangled qubits in the = (||) |t)— |t) [|))/\/2 Bell state. Thus from N 
initial qubits Alice and Bob can prepare N/2 Bell states. Alice sends her A^/2 qubits of the entangled pairs through 
the Pauli channel, which transforms the entangled state into the mixed state 

p' = c(iV'-)(v^-|) = Aiir)(n + A2i0+)(0+i + A3iv^+)(V'+i 

+ (1- Ai -A2 -A3)|V")(V'-|, 

where = (| i) | i) ± | t) | T)) /V2 and = (| i)| T) ± | T)| i)) /V2 are the four BeU states [|. Bob performs 
A^/2 Bell measurements, i.e., projects each of the A^/2 pairs of qubits onto one of the four Bell states with probabilities 
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P(|0-)) = Ai, P(|<^+)) = A2, P(|^+)) = A3, 
P(|i/.")) = 1-Ai-A2-A3. 
Consequently, the estimated parameter values are directly given by 



'1 xcst ^ ^2 , est ^ ^3 

iV/2' 2 ^/2' ^ iV/2' 



We stress that now the estimates of the parameters A,; are always nonnegative, contrarily to the estimation scheme 
with single qubits. 

The average error c^'^"' of the entangled estimation based on the statistical cost function Cg reads |^ 

{N/2)\ 



il+i2+i3+ii=N/2 
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fc=l 

[Ai(l - Ai) + A2(l - A2) + A3(l - A3)] 



N/2 

As we have shown already in |^ c/"' is always smaller or equal to c/'^f, 

A,(iV, A) = {n, a) - cr' {n, a) 

= ^ [5(1 - Ai - A2 - A3)(Ai + A2 + A3) + A1A2 
+A2A3 + A3A1] > 

If we use the fidelity-based cost function c/ we can also write down the average error Cj"' for the entangled 
estimation scheme. It reads 

^ zi!z2!?3!i4! 

il+i2+t3+U=N/2 

X Xl'X'^'XlHl -X1-X2- X3y^CfiN,X) (25) 
and can be evaluated numerically. In figure ^ we show the difference 

Af (N, A) = c/'^P (^N, a) - c/"* (^N, a) (26) 

between the average error obtained with single qubits and entangled qubits. We compare the two average errors when 
the same number of qubits are used. The figure shows that the use of entangled pairs always leads to an enhanced 
estimation and therefore we can consider entanglement as a nonclassical resource for this application. 
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Before ending this section we want to add some comments about our findings. We have found that the mean 
statistical error has the form 

r i 

c,(7V,A) = -^A,(l-A,) (27) 

i=l 

whenever the estimates of the parameters \i are directly given by the frequencies of measurements. This occurs also 
for the relevant case of the entanglement-based protocol for the Pauli channel, but not for the qubit-based protocol. 
Although the parametrization we have used is better indicated since the represent the probabilities of occurence 
of the logic errors, one might be tempted to use a different parametrization A' that gives again the expression ([2^). 
Indeed, this can be done and actually the errors for the A' are smaller. Nonetheless, one can show that in spite of this 
improvement, the scheme based on entangled pairs still gives an enhanced estimation even with the new parameters 
A'. This is not the case of the one-parameter channels, where the use of entangled pairs does not lead to enhanced 
estimation. 



VI. THE GENERALIZED PAULI CHANNEL 



In the previous section we have proposed an entanglement-based method for estimating the parameters which 
define the Pauli channel. As we shall show in this section, this method can be easily extended to the case of quantum 
channels defined on higher dimensional Hilbert spaces. 

Let us start by considering the most general possible trace preserving transformation of a quantum system described 
initially by a density operator p. This transformation can be written in terms of quantum operations Ai as 

p' ^Cp = Y,KA,pA\ (28) 

i 

with \iAiA\ — 1. According to the Stinespring theorem [|l0| this is the most general form of a completely positive 
linear map. The set of operators Ai can be interpreted as error operators which characterize the action of a given 
quantum channel onto a quantum system and the set of parameters A^ as probabilities for the action of error operators 
A,. 

In particular, we can consider the action of the quantum channel, equation (p8|), onto only one, say the second, of 
the subsystems of a bipartite quantum system. For sake of simplicity we assume that the two subsystems are Z?-level 



systems. If only the second particle is affected by the quantum channel, equation (28), then the final state is given by 



p' ^^X^{l®A,)p{l®Ai)^ (29) 

i 

where 1 is the identity operator. In the special case of an initially pure state, i. e. p — (■i/'l, the final state p' 
becomes 

p' = ^A,|^.)(V'.| (30) 

i 

where we have defined \ipi) = (l(g)Ai) {ip). 

Our aim is the estimation of the parameter values Ai which define the quantum channel. This can be done by 
projecting the state of the composite quantum system onto the set of states As in the previous protocols, the 

probabilities A^ can be inferred from the relative frequencies of each state and the quality of the estimation can 
be quantified with the help of a cost function. However, these states can be perfectly distinguished if and only if they 
are mutually orthogonal. Hence, we impose the condition 

(V'.l Vj) = (^1 imA,)^ (mA,) \^) = 5,,j (31) 

on the initial state 1-0). It means that the initial state IV') of the bipartite quantum system should be chosen in 
such a way that it is mapped onto a set of mutually orthogonal states j'i/'i) by the error operators Ai. It is worth to 
remark the close analogy between this estimation strategy and quantum error correction. A non-degenerate quantum 
error correcting code B corresponds to a Hilbert subspace which is mapped onto mutually orthogonal subspaces 
under the action of the error operators. In this sense, a state j?/;) satisfying the condition ( ^l| ) is a one-dimensional 
non-degenerated error correcting code. 
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A necessary but not sufficient condition for the existence of a state 1-0) satisfying equation (|3l|) is given by the 
inequahty 

n<D' (32) 

where n is the total number of error operators which describe the quantum channel and D' is the dimension of the 
Hilbert space. This inequality shows why the use of an entangled pair enhances the estimation scheme. Although the 
first particle is not affected by the action of the noisy quantum channel, it enlarges the dimension of the total Hilbert 
space in such a way that the condition (|l|) can be fulfilled. It is also clear from this inequality that if the state |V') 
is a separable one, a set of orthogonal states can only be designed in principle if the number of error operators n is 
smaller than D' . For instance, in the case of the Pauli channel for qubits considered in the previous section, condition 
( ^ ) does not hold if we use single qubits for the estimation. The action of the Pauli channel is defined by a set of 
n = 4 error operators (including the identity) acting onto qubits {D' = 2). 

Let us now apply our consideration to the case of a generalized Pauli channel. The action of this channel is given 

by 

D-l 

P'= Kf^U^^ppUlp. (33) 

a,/3=0 

where the D"^ Kraus operators Ua^p are defined by 

C/„,^ = {XofiZoT. (34) 

The operation 

is the generalization of a phase flip in D dimensions. The I?-dimensional bit flip is given by 

= l(j-l) mod D). 

These errors occur with probabilities Xa.p, which are normalized, X)a /jio "^"'/S ~ "^^^ states = 0, ...,D — 1 
form a £)-dimensional basis of the Hilbert space of the quantum system. This extension of the Pauli channel to higher 
dimensional Hilbert spaces has been studied previously in the context of quantum error correction [ pT[p^ , quantum 
cloning machines and entanglement . 

Now we will apply the estimation strategy based on condition ( |3l| ) to this particular channel. Thus, we must look 
for a state of the bipartite quantum system which satisfies the set of conditions (|3l|), for the error operators U^^is, i.e. 

It can be easily shown that such a state is 

1 

1^0^") = ^El*)i® 1^)2 (35) 



2 = 



where and \i)2 , j = 0, D — 1, are basis states for the two particles. In fact, the action of the operators 1 Ua,f3 
onto state |'0o,o) generates the basis 



1 ^"^ 

l^^^p) = 1 ® U^J^Pom) ^^J2 e^^"fc|fc), |fc _ 



fe=0 

12 



of D maximally entangled states for the total Hilbert space. Thereby, the action of the generalized Pauli channel 
yields 

D-l 
a,p=0 

when the initial state is p = |'(/'o,o)(''/'o,o|- The particular values of the coefficients A^,/? can now be obtained by 
projecting onto the states |5'a,/3)- The quality of this estimation strategy according to the statistical error is given by 

1 

Cs{N,X) = — — ^ Aq,/3(1 - Aa,/3)(1 - SafiSp.o) 
Q,/3=0 

with N/2 the number of pairs of quantum systems used in the estimation. 
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VII. CONCLUSIONS 



We have examined the problem of determining the parameters that describe a noisy quantum channel with finite 
resources. We have given simple protocols for the determination of the parameters of several classes of quantum 
channels. These protocols are based on measurements made on qubits that are sent through the channel. We 
have also introduced two cost functions that estimate the quality of the protocols. In the most simple protocols 
measurements are performed on each qubit. We have also shown that more complex schemes based on entangled pairs 
can give a better estimate of the parameters of the Pauli channel. Our investigations stress once more the usefulness 
of entanglement in quantum information. 
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